Ferromagnetic Spin Coupling as the Origin of 0.7 Anomaly in Quantum Point Contacts 
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We study one-dimensional itinerant electron models with ferromagnetic coupling to investigate the 
origin of 0.7 anomaly in quantum point contacts. Linear conductance calculations from the quan- 
tum Monte Carlo technique for spin interactions of different spatial range sugg est that 0.7(2e^/^) 
anomaly results from a strong interaction of low-density conduction electrons to ferromagnetic fluc- 
tuations formed across the potential barrier. The conductance plateau appears due to the strong 
incoherent scattering at high temperature when the electron traversal time matches the time scale 
of dynamic ferromagnetic excitations. 



Quantum point contacts (QPC) are narrow constric- 
tions inside two-dimensional electron gas. They con- 
struct one of the building blocks of submicrometer de- 
vices such as quantum dots and qubits [l], Q. The dc 
conductance through a QPC is quantized in steps of 
However, experiments also reveal 
the appearance of an additional shoulder in the conduc- 
tance measurement near O.TGo widely referred to as the 
0.7 anomaly 0]. The origin of 0.7 anomaly in QPC has 
remained a puzzle over almost a decade. The evolution 
of the 0.7Go plateau to O.SGq with magnetic field and the 
enhancement of the g factor [2] have strongly suggested 
that the origin of the anomaly is the electron spin. 

A number of scenarios have been proposed, such as spin 
polarization of the itinerant electrons jh, 16;, 7,, 8 , 9 , iiol 11 1 , 



ferromagnetic correlation |12l . Il3l . Il4| . formation of a 
spin 1^2 magnetic moment in the conductance chan- 




and Kondo effect 



3, [13, 113, Hub- 



lization and antiferromagnetism 27|. These approaches 
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have produced a wide range of different phenomenolo- 
gies, sometimes inconsistent with experiments, and there 
is no widely accepted microscopic theory to date. The 
problem is partly due to the approximate methods used 
in the strongly interacting limit and therefore it becomes 
essential to perform exact calculations to test microscopic 
models against experiments. Here we use numerically 
accurate quantum Monte Carlo technique to study the 
strong correlation effects in QPC devices. 

We find that the 0.7 anomaly at high temperature 
arises from the incoherent electron scattering from itin- 
erant ferromagnetic fluctuations near the Stoner insta- 
bility [2^ , in the strong correlation limit of low electron 
density created by spatially inhomogeneous gate poten- 
tial. We show, through a comparison with a model with 
on-site interactions, that the relevant electron scattering 
is due to the spin fluctuations which are spatially coher- 
ent across the potential barrier. With decreasing temper- 
ature, the magnetic excitation becomes slower than the 
itinerant electrons. The current is then carried by the 
quasiparticles and the 0.7 plateau gradually disappears. 
With the Zeeman magnetic fleld, the 0.7 plateau evolves 



Gate voltage ~ Vg/cosh"(x /Lg) 
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Interacting block 



FIG. f : (color online) Profile of the fD chain with the inter- 
acting block (including seven sites in this figure). The gate 
voltage potential V(x) acts as an adiabatic potential barrier 
through the QPC region. A bias with a ramp passing through 
the QPC is applied across the chain to compute the dc con- 
ductance in the linear response regime. 



to a robust 0.5 plateau in agreement with experiments. 

We model our system using a one-dimensional (ID) 
electron gas with spin-spin interaction among itinerant 
electrons as depicted in Fig|T] with the assumption that 
the higher ID subbands play indirect roles for the first 
conductance plateau. The Hamiltonian reads 



+ j dxi^K^{x)[s{x)f + ]^K2{x)[d,s{x)f^ , (1) 

with = [V't(^)j V'iC^;)]"^ the field operator vector, 

/i the chemical potential, H the Zeeman magnetic fleld, 
(7 the Pauli matrices and V{x) the external gate volt- 
age barrier in order to pinch off the electron current 
through the QPC. V{x) is deflned in our model as 
V{x) = Vg / cosYv' {x / L g) with x — Q corresponding to 
the center of the chain, Lg a characteristic length and 
Vg the gate voltage. The operator s[x) = ip\x)^4i{x) 
represents the spin density of itinerant electrons along 
the chain. Spins interact locally with the coupling con- 
stant Ki{x) — a{x)Ki [Ki < for repulsive on-site 
Coulomb interaction) with a{x) = 1/ cosh{x/ Lg) an at- 
tenuating function with characteristic length Lg. We set 
the spin coupling to adiabatically fall off to take into ac- 
count the screening effects in the leads and to reduce the 
backscattering due to interaction away from the QPC 
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constriction. K2{x) = a{x)K2 {K2 > 0) is the coeffi- 
cient of the gradient term accounting for the ferromag- 
netic Heisenberg interaction. We discretize the contin- 
uum model Hamiltonian to a tight-binding chain of lat- 
tice constant Ax with the nearest-neighbor hopping t. 
Defining t = M - M - VA^ ^{x^) = c„ 

•^0 - -(li + i^), Ji^ ^ and s; = 41 Cp the dis- 
cretized Hamiltonian reads 

H = -i ^ c\^Cj„ - XI " ^ + \^^) 

- X (Jo^pSp • Sp -t- Jittpttp+iSp • Sp+i) , (2) 

with H taken along the z direction and a = ±1 the 
spin index. The microscopic parameters K\ and 
are unknown and we treat Jo, Ji > (for ferromagnetic 
coupling) as the model parameters throughout this let- 
ter. The index p runs only within the interacting block 
in Figd] near the QPC saddle point. The discretiza- 
tion is valid since we are in the low-density limit with 
{cl^Cia) < 1 inside the interacting block. Near the pinch- 
ofF gate voltage, (cl^Cia) at the top of the potential bar- 
rier rapidly approached zero in the following calculations. 
Using m « 0.067rne for GaAs and Ax « 20 nm (exper- 
imental QPC length is about 200 nm, roughly lOAx for 
Lg — 4), t — 1.4 meV. Due to CPU limitations, we 
restrict the interacting block to about seven sites. 

In our calculations, we modify Eq.Q by allowing the 
nonlocal part of the interaction term to extend beyond 
nearest neighboring sites. We define the block spin oper- 
ator S = J2p£biock o^p^p and rewrite a new Hamiltonian 

<ij>,(T icr 

~ (-^o - y j ttpSp ■ sp - y5- 5. (3) 

peblock ^ ' 

Compared to Eq.([21), Eq.® incorporates stronger spin 
interaction among all the spins within the interacting 
block. This modification makes the decoupling scheme 
in quantum Monte Carlo more efficient. Interactions 
beyond nearest neighbors can be thought of a coarse- 
grained effective Hamiltonian on the discretized lattice 
in the low wave- vector limit. The effective interaction re- 
sults from virtual fluctuations to high momentum states 
which are excluded in the discretized model and it takes 
a form similar to the RKKY interaction. Since we are in- 
terested in the low-density limit near the pinch-off regime 
with the effective Fermi wave- vector fcp.off inside the con- 
striction approaching zero, the kp^csRi factor for position 
Ri inside the constriction also goes to zero and the effec- 
tive spin interaction over the interacting block becomes 
predominantly ferromagnetic. 

We use a continuous Hubbard-Stratonovich decoupling 
for the S ■ S term in Eq.([3]) and discrete decoupling for 
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FIG. 2: (color online) (a) dc conductance as a function of 
gate voltage Vg at different values for Lg, Ls, Jo and Ji at 
T — 0.025f. Conductance plateaus form at 0.6 — 0.7 times 
of Go ~ 2e^ /h. (b) The evolution of the plateau to 0.5Go 
with the Zeeman magnetic field, (c) The gradual increase of 
conductance with decreasing temperature, (d) Purely local 
model (Ji = 0) as a function of temperature and the Zee- 
man magnetic field. Local interaction produces qualitatively 
different transport mechanism from the nonlocal model. 

local Sm ■ Sm term. We calculate the dc conductance 
using the Kubo formula in the linear response regime, 

Gdc{^ = 0) = lim^^oRe z e'-*'( [j(t'), i?sd] , 

where ]{t') = lei Ej4,.(i')co,.(i') - cS_,(t')ci,.(i')] is 
the current operator evaluated at the center of the QPC 
and Hsd — — e „ ^{^m)nm,a is the external perturba- 
tion across the chain with V{xm) the normalized source- 
drain bias with maximum (minimum) voltage 1/2 (—1/2) 
on the left (right) hand side as depicted in Fig[T] The 
range of summation in i?sd, \'m\ ^ 100, produced well- 
converged conductance. Conductance by the Kubo for- 
mula is obtained in terms of the bosonic Matsubara fre- 
quencies {ivn = where /? = l/fc^T with the 
Boltzmann constant and positive integer n) and it needs 
to be analytically continued to real frequency to take the 
dc limit Gdc{^ = 0). This task is done by fitting the con- 
ductance defined on the Matsubara frequency into the 
Lehman representation G(w„) — i duj' f'^^J^^, with 
the spectral function p{uj) as the fitting parameter. This 
method has been extensively tested to an excellent agree- 
ment in comparison with the rational function fit [26|. 
After taking the analytic continuation iVn —^uj + ii], we 
obtain the conductance Gdc(w = 0) = p{0)- 

Figllja) plots the dc conductance as a function of 
the gate voltage Vg for different values of Lg and Lg 
at several Jq and Ji values and fixed chemical poten- 
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••T = 0.0141 

T = 0.0251 

T = 0.0331 
AiT = 0.0141, 1^ = 3^=0 

T = 0.0251, J = J = 



II = 0.4t, L = 4. L = 2.5 
J„ = 0.44t, J =0.441 




L,^ ^ ^ ^ ^ ^ ^ ^ . .y ^ l^ = ^-f-^^^-=^- . 
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FIG. 3: (color online) (a) Static ferromagnetic spin correla- 
tion function versus the gate voltage at different tempera- 
tures. The susceptibility has a strong many-body enhance- 
ment (compared to the noninteracting model), (b) Spectral 
function for dynamic ferromagnetic spin susceptibility pFM(t^) 
at different values of gate voltage at T = 0.025t. The dashed 
line indicates the excitation frequency. 



tial /i — OAt when H = Q. With increasing interac- 
tion strength, the plateau evolves from near O.SGo f25j 
to higher values 0.6 — O.TGq in the strongly interacting 
limit with Jo, Ji > /i. Figl^fb) exhibits the gradual evo- 
lution of the plateau near 0.7Go to 0.5Go as the Zee- 
man magnetic field is applied. The significantly wide 
gate voltage interval for the plateau region (AVg w OAt 
at 0.5Go for H = Q.2t) compared to the bare Zeeman 
splitting (AVg = H ^ 0.2t) is a clear indication for the 
enhanced g factor as seen in the experiment Q . FiglSJ^c) 
shows that, with decreasing temperature, the conduc- 
tance plateau consistently moves upward with decreas- 
ing width. Despite the phase problem of the quantum 
Monte Carlo method for temperature T < 0.014t, we are 
able to capture the correct trend as seen in the experi- 
ment for O.OUt <T< 0.033t. T = 0.025t corresponds 
to T « 0.41K, falling within the experimental range [l^ . 
Features in Figl2ta)-(c) correctly reproduce the exper- 
imental results on the temperature and magnetic field 
dependence. 

FigHUd) plots the conductance in the purely local limit 
(Ji — 0). The local spin model is equivalent to the re- 
pulsive Hubbard model 2^, 21, 2^, 2^ with the on-site 
Coulomb parameter U given as J7 = 3Jo/4 by redefining 
the gate potential Vi to absorb the one-body terms. Al- 
though the local limit produces a well-defined 0.7 feature 
at zero field, it is inconsistent with the 0.7 phenomenol- 
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FIG. 4: (color online) (a) pfm{u}) at Vg = 0.5t for different 
temperatures. The peak frequency cUpeak corresponds to fer- 
romagnetic excitation energy, (b) The excitation energy as a 
function of gate voltage Vg at different temperatures. When 
the traversal time rtr (dashed line) is long inside the constric- 
tion, i.e. (jjpeakTtr > 1, the elcctrou scattering becomes strong 
and the conductance plateau results. 



ogy First, the 0.7 feature becomes more pronounced 
as temperature is lowered. Second, more interestingly, 
at finite magnetic field two plateaus appear with the 
0.7 feature shifted to higher conductance and another 
plateau emerging near G ^ 0.3Go. It is very interest- 
ing that these results are consistent with the scenario of 
spin singlet-triplet formation discussed in Refs. 18, |2^ . 
The main difference here is that the spin is self-generated 
from itinerant electrons in our model, not as an exter- 
nal spin [isj or from a quasibound state (2^. Near the 
pinch-off, electron and spin densities are low and the spin- 
singlet does not form and the conductance plateau does 
not appear at H = 0. However, at finite H, the spin 
moment becomes enhanced enough to produce the spin- 
singlet conductance plateau. 

In the presence of nonlocal interactions, an itinerant 
electron interacts with many neighboring sites and the 
resulting spin multiplets are not necessarily S" = or 
5 = 1. With the finite interacting range, a spatially 
coherent ferromagnetic state extends over the interacting 
block at finite H and becomes harder to be flipped by 
electron scattering. Therefore, the nonlocal interaction 
blocks the minority-spin band and the 0.5 plateau results 
from spin splitting, instead of the 0.3 plateau through 
the singlet formation. It has been shown previously that 
ferromagnetic coupling beyond local interaction stabilizes 
the ferromagnetic phase in a uniform ID chain [l3 |. 

Due to the low electron density within the constriction 
at the pinch-off, the ferromagnetic spin correlations are 
considerably enhanced at low T. We compute the ferro- 



magnetic spin correlation function (FSCF), C^p^]^,j ' (vn) = 

W^j:p,ramockIo^"'"^{S^i^)S;nmdr, with l^n = '-f^ 

{n > 0) the bosonic Matsubara frequency, r the imag- 
inary time and N the number of sites inside the inter- 
acting block in Fig[T] FiglHJa) plots the static FSCF 
{I'Q = 0) as a function of Vg at different T values corre- 
sponding to the conductance curves in Figl^l^c). In com- 
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parison with the noninteracting system, the static FSCF 
is significantly enhanced as T is lowered when there is 
interaction. Enhancement of the static FSCF is effec- 
tive only around the pinch-off due to the singular na- 
ture of the ID density of states. The gate voltage range 
for the maximum enhancement of static FSCF also co- 
incides with that of the plateaus in Fig[2ja)-(c), indi- 
cating the effect of strong ferromagnetic correlations on 
the appearance of the plateau. The dynamic ferromag- 
netic spin susceptibility can be obtained by analytically 
continuing the FSCF using the same method employed 
for the conductance. Through the fluctuation-dissipation 
theorem, pfm{'^)i the spectral function for magnetic ex- 
citations can be obtained from the analytically contin- 
ued C^pM^'iivn to + ir]) from positive i^n- In Fig[3fb), 
Pfm{^) has been plotted for different values of Vg at 
T = Ombt corresponding to Figllta). The height of the 
peak grows as Vg increases towards the pinch off, reaching 



its maximum for V„ « 0.4t— 0.5t, around t 



^peak 



0.1135i 



as indicated by the dashed line. As Vg is further in- 
creased, the height of the peak decreases in agreement to 
the decrease in the static FSCF in Fig.[3Ua). However, the 
location of the peak continues to stay at Upeak ~ 0.1135i 
as the characteristic excitation energy at T = 0.025 up to 
Vg « 0.7t where ferromagnetic correlations have almost 
been obliterated as shown in FigOJa). 



Tokura and Khaetskii [ifj addressed the effect of scat- 
tering for electrons in the current due to the local param- 
agnons as the origin of the 0.7 anomaly in the QPC using 
the second order perturbation approach with local inter- 
actions. Following Biittiker and Landauer in Ref.js^, we 
argue that magnons, perceived as spin waves with charac- 
teristic frequency ujpeak , can strongly interfere with elec- 
trons while transmitting through the gate voltage barrier 
at frequencies Upeak ~ ^/Ttr where Ttr is the traversal 
time for tunneling through the barrier in the absence of 
magnons defined as J^^ ^^^2[F(" ™)-_e] ' ^2 
the classical turning points, E fi the incident energy 
and V{x) = Vg / cosYi' {x / L g) the gate voltage barrier. 

Figlllja) presents pfm{'^) for Vg — Q.ht (maximum 
FSCF in FigEl^a)) at three different T values. The char- 
acteristic frequency ujpeak decreases as the temperature is 
lowered and therefore the spin waves soften at the onset 
of Stoner instability in the system as explained in detail 
in Ref. [1^ . FigUJ^b) plots ujpeak as function of Vg for these 
three different temperature values along with the recip- 
rocal traversal time 1 / Ttr (nearly constant as a function 
of Vg). While LOpeak ^ ^/Ttr^ strong interference between 
the spin fluctuations and electrons tunneling through the 
gate voltage barrier leads to the suppression of the con- 
ductance for the gate voltage values in the vicinity of the 
pinch off. The manifestation of this suppression is the 
appearance of a plateau in the conductance as a function 
of the gate voltage near the pinch off. With decreas- 
ing temperature and therefore LOpeak ^ l/''"tr, as seen in 



FiglUa) and (b) for T = 0.014^, the itinerant electrons 
begin to see an effectively static ferromagnetic mean field 
and scattering of electrons off magnetic excitations be- 
comes coherent. As a result, the conductance increases 
and the plateau disappears. 

We have demonstrated the incoherent electron scatter- 
ing due to ferromagnetic spin-spin correlations in a quasi- 
ID chain as the underlying physics behind 0.7 anomaly 
phenomenon in the QPC. We thank Jonathan Bird and 
Igor Zutic for helpful discussions. This project was sup- 
ported by NSF DMR-0426826 and we acknowledge the 
CCR at the SUNY Buffalo for computational resources. 
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